Equivalent forms of Dirac equations in curved 
spacetimes and generalized de Broglie relations 

^ ; Mayeul Arminjon^ and Frank Reifler^ 

^ . ^ Laboratory ''Soils, Solids, Structures, Risks" 

1—5 ■ (CNRS and Universites de Grenoble: UJF, G-INP), 

■ BP 53, F-38041 Grenoble cedex 9, France. 
^\ 2 Lockheed Martin Corporation, MS2 137-205, 

199 Borton Landing Road, Moorestown, New Jersey 08057, USA. 

' Abstract 
> '■ 

. One may ask whether the relations between energy and frequency and 

I between momentum and wave vector, introduced for matter waves 

' by de Broghe, are rigorously valid in the presence of gravity. In this 

. paper, we show this to be true for Dirac equations in a background 

[ of gravitational and electromagnetic fields. We first transform any 

^— I ' Dirac equation into an equivalent canonical form, sometimes used in 

■ particular cases to solve Dirac equations in a curved spacetime. This 
! canonical form is needed to apply the Whitham Lagrangian method. 

^ I The latter method, unlike the WKB method, places no restriction on 

H ■ the magnitude of Planck's constant to obtain wave packets, and fur- 

- - - thermore preserves the symmetries of the Dirac Lagrangian. We show 

using canonical Dirac fields in a curved spacetime, that the probabil- 
ity current has a Gordon decomposition into a convection current and 
a spin current, and that the spin current vanishes in the Whitham 
approximation, which explains the negligible effect of spin on wave 
packet solutions, independent of the size of Planck's constant. We fur- 
ther discuss the classical-quantum correspondence in a curved space- 
time based on both Lagrangian and Hamiltonian formulations of the 
Whitham equations. We show that the generalized de Broglie rela- 
tions in a curved spacetime are a direct consequence of Whitham's 
Lagrangian method, and not just a physical hypothesis as introduced 
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by Einstein and de Broglie, and by many quantum mechanics text- 
books. 

PACS numbers: 03.65.Pm Relativistic wave equations 
03.75. - b Matter waves 

04.62. + V Quantum fields in curved spacetime 

1 Introduction 

Some might argue that quantum mechanics became a universal theory of mat- 
ter when, at the last Einstein-Bohr debate at the Solvay Conference in 1930, 
Einstein proposed a weight measurement to observe unobtrusively a particle 
decay in order to contradict the Heisenberg energy-time uncertainty relation 
[1] . After nearly being defeated by Einstein in the debate, Bohr surprisingly 
countered with a general relativistic gravitational argument. Henceforth, the 
relation between gravity and quantum mechanics was to become an impor- 
tant question in fundamental physics. 

It could also have been questioned whether the relations between energy 
and frequency and between momentum and wave vector, introduced for mat- 
ter waves six years earlier by de Broglie [1] , were rigorously valid in a general 
curved spacetime. This question can be shown to be equivalent to the ques- 
tion of whether sufficiently small wave packets travel along classical paths 
consistent with the de Broglie relations. Recall that a wave packet is a wave 
whose amplitude, frequency, and wave vector vary slowly over a region of 
spacetime comparable to a period or wave length. (E.g., for an electron trav- 
eling at half the speed of light, the wave length is approximately 5 x 10~^^ 
m, and a typical wave packet has dimensions 10"^ m [2].) The de Broglie 
relations can be observed for such wave packets. 

The well known WKB approximation is commonly used to derive wave 
packet approximations in quantum mechanics (e.g., [3]). The WKB approx- 
imation is based on taking the limit as a physical constant, namely Planck's 
constant h, approaches zero. However, the spin connection in the Dirac 

^ More generally, a common conception is that classical physics emerges from quantum 
mechanics in the limit as Planck's constant h approaches zero. However, the limit ft — > 
"is not well defined mathematically unless one specifies what quantities are to be held 
constant during the limiting process" [4]. It is interesting to note that the most classical 
behaving Gaussian wave functions, the coherent states of the ordinary harmonic oscillator, 
whose expected position and momentum obey classical equations by Ehrenfest's theorem, 
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equation of a curved spacetime has no effect in the first WKB approximation 
(i.e., the one retaining only the zero order term in h) [5]. The assumption 
that the spin connection can be neglected, as it would be in a first WKB 
approximation, is unnecessary and is too strong for many applications in a 
curved spacetime, or even in a Minkowski spacetime with arbitrary coordi- 
nates. Note that throughout this paper, except for the brief description of a 
post Newtonian approximation in Section [31 we may set both the speed of 
light c and Planck's constant h equal to one. 

The Whitham approximation [9], which we adopt in this paper, places 
no restriction on Planck's constant. To implement the Whitham approxima- 
tion and to show that it leads to propagation along classical paths, we will 
first show in Section [2] that any Dirac equation in a curved spacetime can 
be transformed into an equivalent canonical form known in the literature as 
the "local representation" |6], [7]. In general, transformation to equivalent 
canonical form is a necessary step to simplify a Dirac equation so that prop- 
agation along classical paths can be derived. It will be evident in Theorem 1 
of Section [2] that Planck's constant h does not appear in the transformations 
mapping Dirac equations to their equivalent canonical forms [Hj, [7]. Pre- 
viously, these canonical forms, or "local representations" as they are called 
in the literature, have only been discussed in the special case of orthogonal 
coordinates [6], [7], [8]. 

Then in Section [31 with each Dirac equation transformed into equiva- 
lent canonical form, we apply Whitham's Lagrangian method [5] to derive 
wave packets in general curved spacetimes. Whitham's method preserves the 
symmetries of the Lagrangian, and in particular, the Whitham wave packet 
equations conserve the probability current. We also show that generalized de 
Broglie relations, as well as COW and Sagnac type terms [TU], emerge from 
the Whitham equations. It will become clear in Sections 2 and 3 that for 
every Dirac equation transformed into equivalent canonical form, the general- 
ized de Broglie relations have no other meaning than the fact that sufficiently 
small wave packets propagate along classical paths in a background of grav- 
itational and electromagnetic fields. This is what is observed in experiments 
and therefore more physically precise than the statement often made that 

do not resemble wave packets in the limit fi, — >■ [5] . 

Note, however, that the electron's magnetic moment, predicted by the Dirac equation, 
is not contained in the wave packet approximation. Indeed, quoting from Ref. "The 
uncertainty principle, together with the Lorentz force, prevents spin-up and spin-down 
electrons from being separated by a macroscopic field of the Stern-Gerlach type." In 
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particles with a given energy and momentum possess a frequency and wave 
vector given by the de Broghe relations. In fact, the generalized de Broglie 
relations are a direct consequence of Whitham's method applied to each Dirac 
equation transformed into equivalent canonical form, and not just a physical 
hypothesis as introduced by Einstein and de Broglie, and by many quantum 
mechanics textbooks. 

In the WKB approximation of the standard Dirac equation in a curved 
spacetime, classical trajectories are derived from the Gordon decomposition 
of the probability current = + into a convection current and 
a spin current J^f [3] . In Section [3l we also prove the existence of the Gor- 
don decomposition for Dirac equations transformed into canonical form. We 
further show that in the Whitham approximation, the spin current van- 
ishes, which explains the negligible effect of spin on wave packet solutions, 
independent of the size of Planck's constant h. 

It is also clear that the canonical forms (or "local representations" ) of the 
Dirac equations, while not unique, are the preferred representations to un- 
derstand certain phenomena associated with the Dirac equations in a curved 
spacetime, particularly, the emergence of classical physics and its conserva- 
tion laws in a quantum world. Section |4] concludes this paper with a discus- 
sion of the classical-quantum correspondence in a curved spacetime based on 
both Lagrangian and Hamiltonian formulations of the Whitham equations. 
In this section we also include results from a previous analysis of the classical- 
quantum correspondence [T2j, which can be applied to the canonical forms 
of Dirac equations in a curved spacetime considered in this paper. 

practice, wave packet splitting in Stcrn-Gerlach experiments is only observed using neutral 
atoms or molecules, wliieli are undisturbed by the Lorentz force [3]. In Seetion[3J the wave 
packet approximation is expressed by neglecting in the Lagrangian the variation in the 
amplitude of the wave function as compared to the variation of its phase. This leads to 
wave packet equations which do not involve the electron's magnetic moment. 
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2 New Representations of the Dirac Equa- 
tion in a Curved Spacetime and their Equiv- 
alent Canonical Forms 

Shortly after Dirac discovered his celebrated four component wave equation: 

rd,^> = (1) 

together with its conserved probability current: 

= c^+A7^^, (2) 

his equation was studied in its widest representations for a Minkowski space- 
time [13], [H]. 

In Eq. ([1]), the Dirac field is a four component complex function of 
spacetime coordinates x^^, fi = 0,1, 2, 3, whose partial derivatives with respect 
to are denoted as The Dirac gamma matrices 7^, acting on satisfy 
the anticommutation formula: 

+ YY = '2v^"U, (3) 

where rj^^'^ is the inverse of the Minkowski metric tensor rj^^, and I4 denotes 
the identity matrix acting on the Dirac field The mass, the speed of light, 
and Planck's constant are denoted by m, c, and h, respectively. Repeated 
indices are summed. 

In Eq. ([2]), '^'^ denotes the complex conjugate transpose (Hermitian 
conjugate) of the Dirac field \Ef, and A is a hermitizing matrix for the Dirac 
gamma matrices 7^^. That is, [I3], [11]: 

A+ = A, 

(4) 

7M+ = A'yi'A-^, 

where 7^*^ and A~^ denote the Hermitian conjugates of the matrices 7^ and 
A, respectively. The hermitizing matrix A is uniquely determined by the 
matrices 7^^ up to a nonzero real scalar multiple |13j . 

For a Minkowski spacetime, assuming that (7^, A) are chosen to be con- 
stant matrices satisfying Eqs. ([3]) and (jl]), every solution of Eq. ([1]) satisfies 
the Klein-Gordon equation, and the probability current in Eq. is then 
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also conserved. That is, Eqs. and (jl]) are the only conditions that the 
constant matrices (7^, A) need satisfy. 

The "coefficient matrices" (7^^, A) in Eqs. ([1]) and ([2]) are far from unique. 
Given a Dirac field and any set of constant coefficient matrices (7^^, A) sat- 
isfying Eqs. ([3]) and (jl]), they may be transformed by any constant complex 
4x4 matrix 5* as follows: 

$ = ^-1*, 

Y = S-'-f^'S, (5) 
A = S+AS. 

Such a transformation 5* is called a "similarity transformation" or a "spin- 
base transformation" , the latter referring to simply a change of basis for the 
four components of the Dirac field It is straightforward to see that Eqs. ([T]) 
— (jlj) are invariant under all similarity transformations S by Eq. Thus, 
in the widest sense no restriction to a smaller group of transformations was 
deemed necessary in the early development of the Dirac equation [13], [H]. 
Nevertheless, further choices were necessary when the Poincare group of co- 
ordinate transformations of the Dirac equation was considered. At least three 
possibilities have been considered for Poincare coordinate transformations in 
a Minkowski spacetime as follows: 

• A) The Dirac field \1/ transforms as a quadruplet of complex scalars 
under a coordinate transformation. [15], [IS], [IZ], [IB] 

• B) The Dirac field \l/ transforms as a complex four-vector under a 
coordinate transformation. [I9] 

• C) The Dirac field \E' transforms as a quadruplet of complex scalars 
under a coordinate transformation, which is then followed by a similar- 
ity transformation. (The combined transformation leaves the constant 
coefficient matrices (7^^,^) invariant.) [T3] 

For a Minkowski spacetime with Poincare coordinate transformations, all 
three possibilities may be considered. However, with general coordinate 
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transformations, as required for a curved spacetime, only the first two pos- 
sibilities (A) and (B) exist. Since in a general spacetime, the possibility (C) 
does not exist, it was replaced by the possibility (A) in what has become the 
standard Dirac equation, which was proposed independently by Weyl [15] 
and by Fock [16], hereafter called the "Dirac- Fock- Weyl" (DFW) equation. 
See Refs. [T7] and [IB]. Possibility (B) was investigated recently [IH], which 
became the impetus for a more general study of Dirac equations [20], [21], 
[22]. 

General Dirac fields of type (A) will be said to belong to the Quadruplet 
Representation of the Dirac theory (or QRD theory). General Dirac fields 
of type (B) will be said to belong to the Tensor Representation of the Dirac 
theory (or TRD theory). It was recently shown that in an open neighborhood 
of each spacetime point, every TRD equation is in fact equivalent to a QRD 
equation and vice versa [20]. Since TRD equations are locally equivalent to 
QRD equations, we will only consider QRD equations in this paper. Note 
that there are QRD equations in a curved spacetime that are not locally 
equivalent to any DFW equation [20] . 

In a further evolution of the Dirac equation, which applies also to a 
Minkowski spacetime, the partial derivatives in the Dirac equation ([1]) were 
replaced by covariant derivatives = 9^ + where F^ are four 4x4 
complex matrices [f] , called "spin connection matrices" , acting on the four 
components of the Dirac field At the same time, the Dirac equation was 
generalized by substituting a general metric Qfj^^, of Lorentz signature and 
determinant g for the Minkowski metric 77^,^ in the anticommutation formula 
03]) for the Dirac gamma matrices: 

+ YY = 2^^^i4- (6) 

For the results in this paper, mild restrictions must be placed on the metric 
components g^u, namely that 5^00 > and the components gjk for j, k = 1,2,3 
form a negative definite 3x3 matrix. Even though these mild conditions 
hold for almost all spacetime metrics gf^i, of interest, the Godel metric is a 
notable exception [23], [21]. 

For this generalization, the coefficient matrices (7^^, A) defining the Dirac 
equation are augmented to become "coefficient fields" (7^, A, F^), which now 

^ In the case of a Majorana representation of the Dirac field, the coefficient matrices 
(7^^, A) are pure imaginary, and the Dirac equation is real. In this case the spin connection 
matrices are real. 
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may vary with the spacetime point. Then, in order that the Dirac equation 
([T]) be invariant, the transformation equations are augmented with the 
transformation of the spin connection matrices under "local similarity 
transformations" S (i.e., similarity transformations 5* that may also vary 
with the spacetime point ) as follows: [T8] 



(7) 

= {d^ + r^) s. 

Indeed, for this type of transformation, we have the covariant derivatives 
transforming as: = o Z}^ o S. Transformations of the kind given by 
Eqs. (E]) and ([7]) we will call local similarity transformations "of the first 
kind" . 

Local similarity transformations "of the second kind" are defined by set- 
ting: 

= r,. (8) 

For this second kind of transformation, given by Eqs. (E]) and (|H]), we have 
the covariant derivatives transforming as: = D^. 

Two Dirac equations will be said to be "equivalent" or "classically equiv- 
alent" if there exists a local similarity transformation \1/ — S'"^\l/ of any kind 
that transforms the solutions of one Dirac equation into the solutions of the 
other. From Eqs. and (jSj), the conserved probability currents for two 
equivalent Dirac equations are equal. Hence in any spacetime, scattering ex- 
periments will give the same results regardless of the representation of a given 
Dirac equation. However, in a first quantized theory, which is our concern 
in this paper, a local similarity transformation S may not intertwine with 
the quantum mechanical operators. In that case, the operators correspond- 
ing to a given observable generally will not have the same eigenvalues before 
and after the transformation S. This has been shown for the Hamiltonian 
(or energy) operator in previous work [25]. Thus, two Dirac equations that 
are equivalent as partial differential equations via a local similarity transfor- 
mation S, need not be equivalent with respect to all quantum mechanical 



Because S depends on the spacetime point, Schliitcrt, Wietscliorke, and Greiner call 
S a "local transformation" [6], [7]. 

^ The notion of equivalence here is somewhat different than the notion used in Ref. 
[25] . where equivalence was applied only to classify the coefficient fields (7'^, A) without 
requiring the existence of a map — >■ S~^"^ between the solutions of two Dirac equations. 
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operators [25], [2B], [27]. 



The standard "Dirac Lagrangian" applies to the DFW equation [T7] and 
has to be extended to include the coefficient field A in the exact place of the 
constant hermitizing matrix valid for DFW [25] : 

(9) 

Note that the Lagrangian (Q is the local expression of a global Dirac La- 
grangian based on a general Dirac operator ^ acting on the cross-sections of 
a spinor bundle E over the spacetime. |^ Thus, the Lagrangian (|9]) gives a 
generalized formulation of the Dirac theory for a general Dirac operator 'p on 
a curved spacetime. We will see in Theorem 1 that the equivalent canonical 
forms of DFW equations require such a generalization. 

The Euler-Lagrange equations for the Lagrangian iQ give the following 
generalized Dirac equation [20], [21] : 



where 



1 ITfl C 

YD,^ + -A~'D, ( AY) ^ = -- (11) 

D^r = v^Y + '^.Y-r'^f. (12) 

D,A = d,A - r+A - AT,, 

and we define the Levi-Civita covariant derivatives acting on the Dirac 
field and the coefficient fields (7^, A) as follows: 



^ A globally defined generalized Dirac Lagrangian has the form: 

^[{i;,piP)~iP^P,ij)+2im{i',i;)]. (10) 

where ^ is a Dirac operator acting on the cross-sections ip oi a spinor bundle E over the 
spacetime, and ( , ) denotes a hermitizing metric on the fibers of E. See Ref. [20], Sect. 
2.1 and the references therein for the precise definitions. Once any coordinate chart of the 
spacetime and any frame field on the spinor bundle have been chosen, one gets the local 
expression of the global generalized Lagrangian as Eq. (|9]). In particular, in Eq. ([9]) and 
in the rest of this paper, is the column vector made with the components of "tp in the 
chosen frame field on E. See Ref. [20l, Sect. 2.2. 



9 



= d,Y+{'p,]r (13) 

where {p^} are the Christoffel symbols belonging to the Levi-Civita connec- 
tion. The covariant derivatives extend to by the formula D^^^ = 
(Dfj^^)~^, and similarly V^^'"'" = (V^\l/)^. As usual, covariant derivatives 
extend to products of fields via Leibniz's rule for differentiating products. 

If Dfj^ (^7^) = 0, then the generalized Dirac equation ( 1TT]1 reduces to 
normal form: 

YD,^ = (14) 

Normal Dirac equations generalize the DFW property that the coefficient 
fields (7^^, A) be covariantly constant: D^'-f'^ = and D^A = 0. One can show 
that the weaker normality condition (^7^) = is preserved under all local 
similarity transformations of the first kind [20]. The normality condition is 
also preserved under all coordinate transformations. 

Thus we have several invariant classes of Dirac equations. First, we have 
the class of Dirac equations for which the coefficient fields (7^^, A) are covari- 
antly constant — that is, D^'-f'^ = and D^A = 0. This class contains the 
DFW equations as a proper subset. Second, we have the class of Dirac equa- 
tions for which (^7^) = 0. This is the class of normal Dirac equations 
f lT^ which contains the first class as a proper subset. Finally, we have the 
class of generalized Dirac equations flTT]) which contains the other two classes. 
Each of these classes is invariant under all coordinate transformations and 
also under all local similarity transformations of the first kind |20j . 

A non-invariant class of Dirac equations that we will discuss in this paper 
is the class of the QRD-0 equations, in which the contracted spin connection 
matrix F = Y'^^, = 0. See Ref. [20], Sect. 3.2.1. Equations in the QRD-0 
class (with F = 0) are generalized Dirac equations which may be written as 
follows: 

rO,^ + \a-W, ( AY) ^ = (15) 

with the Levi-Civita covariant derivatives acting on the coefficient fields 
(7^, A) as previously defined in Eq. If (^7'^) = 0, then the QRD-0 
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equation fllSI) reduces to the simpler normal form: 

rd,^ = (16) 

Note that a QRD-0 equation which is both normal as in Eq. (fT6l) and 
equivalent to a DFW equation is called a "local representation" of the DFW 
equation by other authors [6], [7], [8]. Finding a "local representation" for a 
DFW equation often simplifies the process of deriving its solutions. Previ- 
ously, this "local representation" of the DFW equation in a curved spacetime 
in the form of Eq. f|T6|) has only been discussed in the special case of orthog- 
onal coordinates [6], [7], [Sj. 

Given any normal Dirac equation f lT^ . one can construct an equivalent 
normal QRD-0 equation flTB]) in terms of a basis of solutions of the massless 
equation associated with flTll) : 

7^D^^ = r [d, + r^) ^ = 0. (17) 

This is useful when solutions to the massless equation flTTI) are known, for 
which there are many examples in general relativity, including all diagonal 
metrics and Godel type metrics [28] — [32]. Indeed, consider a local similarity 
transformation 5* which takes the Dirac field to and the coefficient fields 
(7'',A,r^) to {^^'.AJ^, according to Eqs. ([5]) and 0. Recall that the 

contracted spin connection matrix F = 7^r^ = if the transformed equation 
is a QRD-0 equation. Then from Eqs. ^ and ([7]), we get: 

Y{d, + V,)S = ^, (18) 

whereby any four linearly independent solutions to the massless equation 
(fT7|) can be used to form the columns of the matrix-valued field S. In this 
case, an equivalent normal QRD-0 equation can be explicitly and globally 
constructed. 

We can show that any generalized Dirac equation ffTTl) . with very minor 
conditions imposed on the spacetime metric (7^1,, is equivalent to a normal 
QRD-0 equation fITBl) . in an open neighborhood of each spacetime point, by 
applying local similarity transformations of the first and second kind to the 
Dirac field Previously, this so called "local representation" of the DFW 
equation in a curved spacetime has only been discussed in the special case 



11 



of orthogonal coordinates [B], [I], [B]- Here we generalize it to essentially all 
Dirac equations: 

THEOREM 1. Let U be any open subset of a spacetime on which lo- 
cal coordinates are defined. Suppose that the metric components g^p for 
/i, z/ = 0, 1, 2, 3 in U satisfy Qqq > and the components gjk for j, k = 1,2,3 
form a negative definite 3x3 matrix. Then, for any choice of smooth coeffi- 
cient fields (7^^, A) and any choice of covariant derivatives D^^ = acting 
on smooth Dirac fields \l/ defined on U, there exists a smooth local similarity 
transformation 5* of the first kind, composed with a smooth local similarity 
transformation T of the second kind, taking \1/ — > (T o 5*)^^ which trans- 
forms the generalized Dirac equation with smooth coefficient fields (7^, A, T^) 
into an equivalent normal QRD-0 equation, in an open neighborhood of each 
point Xq G U. 

Thus, we may regard the normal QRD-0 equations (fT6|) as canonical forms for 
all generalized Dirac equations f ITT]) . in open neighborhoods of each spacetime 
point. 

The proof of Theorem 1 relies heavily on the theory of linear hyperbolic 
partial differential equations [20], [33], [31]. Since it is not particularly en- 
lightening beyond the explicit construction given above for the normal case, 
we will postpone writing out the full proof of Theorem 1 until the Appendix. 

Note that every DFW equation has the normal form fll4p . By Theorem 
1, every DFW equation is equivalent to a normal QRD-0 equation (called a 
"local representation" of the DFW equation by Schliiter, Wietschorke, and 
Greiner) which generally is not a DFW equation [6], [7], [Sj. Conversely, we 
can show that not every normal QRD-0 equation is equivalent to a DFW 
equation, so that DFW equations are in fact equivalent to only a proper 
subset of the possible normal QRD-0 equations. 

EXAMPLE. Consider the flat metric gf^j, on whose line element ds in 
rotating cylindrical coordinates {t, r, 0, z) has the form: 

ds'^ = (c dtf - dr^ - {d(f) + u dtf - dz^ . (19) 

Using the orthonormal tetrad [35] as indicated by the parsing of the metric 
in Eq. (fT9|) . the DFW equation is given by: 

1 771 C 

YD,^ = Ydf.^ + -7'^ = (20) 
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A simple local similarity transformation of the first kind: 

5"^ : ^ ^' = (21) 

which is independent of Planck's constant, time independent, and also is 
independent of the rotation rate w, transforms the DFW equation (!20|) into 
a Dirac equation of the canonical form (fT6!l . Note that the canonical Dirac 
equation (fT6l) does not have the Mashhoon term [36], [10], whose presence 
or absence thus depends on the chosen representation of the Dirac field. See 
also Ryder [57] . 

3 Whitham's Lagrangian Method — The Main 
Theorem 

In this section we apply Whitham's Lagrangian method [9] to derive wave 
packets for Dirac equations in general curved spacetimes. Whitham's method 
preserves the symmetries of the Lagrangian, and in particular, the Whitham 
wave packet equations conserve the probability current. We show in this 
section that generalized de Broglie relations, as well as COW and Sagnac 
type terms [10] , emerge from the Whitham equations after transforming each 
Dirac equation into an equivalent canonical form. Thus, the normal QRD-0 
representations (or canonical forms) of the Dirac equations are the preferred 
representations to express the generalized de Broglie relations in a curved 
spacetime. 

It is noteworthy that the Whitham approximation places no restriction 
on Planck's constant h. Indeed, the transformation which takes the Dirac 
equation to its canonical form, is independent of Planck's constant h. This is 
obvious in the example above, as seen in Eq. ( 12T|) . In fact, this independence 
is a general fact that follows from Eqs. ( l82l) and ( 184|) used in the proof 
of Theorem 1 in the Appendix. Thus, the Whitham approximation is not 
equivalent to the WKB approximation, since the latter does not require any 
transformation of variables. 

^ As stated in Section [21 two Dirac equations that are classically equivalent, need not 
be equivalent with respect to their quantum mechanical energy- momentum operators |25j , 
|?7] . Clearly, the Whitham approximation also distinguishes them. A striking ex- 
ample is a Dirac equation with a Mashhoon term [T0|, [36], [37]. Applying the Whitham 
approximation directly to a Dirac equation with a Mashhoon term, without first trans- 
forming the Dirac field ^E*, docs not produce wave packet motion along classical paths. 
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Including both gravitational and electromagnetic external fields, the gen- 
eralized Lagrangian for the Dirac equation can be written as follows: 



(22) 

with the covariant derivatives Z)^ defined by: 

D^ = d^ + T^ + -V^, (23) 

where V"^ are electromagnetic gauge potentials and e is the electric charge. 
We can display the interaction terms of the Lagrangian ( l22l) more explicitly 
by expressing the Lagrangian ( l22i) as a sum of a free and an interaction part 
as follows: 

L = ^/^if[\i/+AY (d^^) - (df,^)^ AY^ + ^i<i/+A^] 

(24) 



where = c\I^+747^\I^ is the probability current, and F = 7^r^ is the con- 
tracted spin connection matrix, and also noting that since A is a Hermitizing 
matrix for the Dirac matrices 7'^, we have from Eq. dl]): 

T+A = T^Y'-A = ^+A7^ (25) 

Since we can first transform any Dirac equation into a normal QRD-0 equa- 
tion (or canonical form) as stated in Theorem 1 of Section [21 we can transform 
the fields in this Lagrangian so that (^7'^) = and F = 7^F^ = 0. Thus, 
instead of the Lagrangian (1221) . we may substitute an equivalent canonical 
Dirac Lagrangian: 

L = l^+AY (d^"^) - (d^^)^ AY"^ + ^t^+A^] 

(26) 

Clearly, the normal QRD-0 equation f lT6|) is derived from the Euler-Lagrange 
equations of the Lagrangian fj26|) by setting the electromagnetic gauge po- 
tentials equal to zero and using the normality condition: (^7^) = 0. 

For Whitham's method, we set \Ef = x^*^ where x = X (^) is also a 
complex wave function with four components, and 6 = 6 (X) is a real phase at 
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each point X in the spacetime. Whitham's method assumes that x is slowly 
changing compared to the rapidly changing phase 6, so that we may obtain 
approximate wave packet solutions to the Dirac equations by neglecting d^x 
with respect to {df^6) x in the Lagrangian. Substituting the wave function 
^ = xe*^ into the Lagrangian fl26|) and using this approximation, we get: 



(^-hd,e - -V,^ x^A^X - mcx+Ax] . (27) 

In Whitham's method, this Lagrangian governs the wave packet motion. 
Clearly the Lagrangian f l27|) is invariant under the global gauge symmetry 
6 ^ 6 + T, where r is a real constant. This leads by Noether's theorem to 
the conservation of a current. In Section 13.31 we will derive explicitly the 
conservation of the probability current for this Lagrangian. Thus, our goal 
in this section is to derive the Euler-Lagrange equations for the Lagrangian 
(!27|) . which by change of field variables leads to the following main theorem 
of this paper. 

THEOREM 2. Define a four- vector field and a scalar field J, related to 
the amplitude x ^^'^ phase 9 of the wave function = xe*^ as follows: 

u. = -^dj - ^K, 

= gt^t^uy, (28) 
J = cx'^Ax- 

i) Then the Euler-Lagrange equations for the Whitham Lagrangian (1271) im- 
ply the following equations for the fields and J in a curved spacetime: 

g^.uU^'u" = 1, 
(Ju^^) = 0, 

(29) 

where F^,^ = V^K — VjyV^ is the electromagnetic field tensor. 

ii) The four-vector field is a unit velocity field, such that = Ju^ is a 
conserved probability current. 
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iii) The integral curves x'^ (s) of the four-vector field u^, parameterized by 
arc-length s, are given by the classical equations: 



mc 

along which the scalar field J satisfies: 



If + {/p}«^«'' = -S^^«^ (30) 



ds 

The proof and interpretation of Theorem 2 will occupy the rest of Section 3 



JV^u^. (31) 



3.1 Euler-Lagrange Equations for the Wave Packet La- 
grangian 

The Euler-Lagrange equations for the amplitude x = X (^) phase 6 = 
9 {X) can be derived from the wave packet Lagrangian L in Eq. (I27|) as 
follows. First we have from Eq. (1271) : 

^ = [{-fid.e - ^V,) A^^x - mcAx] , 

(32) 



dL 



-hcyf^x'^ AYx- 



Then, since no derivatives of x"*" and only derivatives of 9 occur in the La- 
grangian (127|) . we set equal to zero, using Eq. ( 132|) . the following expressions: 

= 1^ = ^= [{-nd,e - ^V,) A^X - mcAx] , 



(33) 



SL _ a dL 



() = — = d 



which then gives the following Euler-Lagrange equations: 

{-hd^e - W^) AYX= mcAx, 

(34) 

d^icy/^X^^l'^x) = 0. 

For the wave function = xe*^, with phase 6, the wave covector is K^^ = 9^^. 
Thus u = —Kq = —doO is the angular frequency of the wave. We define a 
four-vector field as follows: 
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p^^-hd^9= -hK^. (35) 

It will be shown in Section |4] that = —P^-, where are canonical momen- 
tum variables. 

Eq. (135!) expresses the generalized de Broglie relations = hK^ between 
the canonical momentum variables P^ and the wave covector K^. We also 
define a four-vector velocity field from the usual classical equation with 
kinetic and potential terms as follows: 

= mcu^ + -V^. (36) 
c 

From Eq. (135|) we have d^py = di^p^. Substituting \1/ = xe*^ into Eq. ([2]) we 
have = cx^^7^X- We denote 7 (m) = m^7^. Then, Eqs. ( 134|) become: 

7(^)X = X, (37) 
d,{V^J'') = 0, (38) 
d^^Pu - d^Pf, = 0. (39) 

The first equation ( 1371) is an algebraic eigenvalue equation. The second equa- 
tion (l38ll can be written as the covariant conservation of the probability cur- 
rent, V^J^ = 0, where denotes the Levi-Civita covariant derivatives. 
Since p^ = —hd^9, the third equation fj39l) expresses the equality of mixed 
partial derivatives of 6. Furthermore, since the left-hand side of the third 
equation ( 139|) is an antisymmetric tensor, the partial derivatives 9^ can be 
replaced by Levi-Civita covariant derivatives V^. Thus, Eqs. (1371)-( 139|) be- 
come the following covariant equations: 

l{u)x = X, (40) 

v^r = 0, (41) 

^^Pu-^uP^. = 0. (42) 

We will show below that these equations, taken together, reduce to a set 
of quasi-linear partial differential equations describing a scalar density J = 
cx'^^X ci-nd the four-vector velocity field u^, whose integral curves are classi- 
cal relativistic trajectories. We will further show that certain of these equa- 
tions give rise to initial conditions, and of the rest, only four equations are 
independent. First, let us derive a dispersion relation from the algebraic 
equation (1101) . 
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3.2 Dispersion Relation 



From the algebraic equation poj) . and the anticommutation relation of Dirac 
gamma matrices (|6]), we have: 

{g^'^u^u^) X = 7 {uf X = 7 (n) X = X- (43) 

Equation fH3|) implies that at any spacetime point where the wave function 
X is not zero, the four- vector velocity field satisfies u^n^ = 1. From Eq. 
this gives the dispersion relation: 

g'"" [p, - -V^) {Pu - -K) - m^c" = 0. (44) 

Since po = e/c, where e is the energy and Pj = —pj for j = 1,2,3 are 
momentum variables, Eq. f H^ is a quadratic equation for the energy e. 

Let us consider the dispersion relation f H^ in the absence of the electro- 
magnetic potentials V^. We have: 

g'^P.P. = g'' (^) ' + 2g'^P, (^) + g^'PjPk = m^c\ (45) 

That is. 



cg^^Pj ± c\J {g^ipjf - 5-°° {g^^PjPk - m^c^) 



(46) 



^00 

Choosing positive energy e > and setting: 

/» = l-^. /' = f, (47) 

where 0, (jy' , and 0-^'^' are gravitational potentials, and 5^^ is the Kronecker 
delta, equal to one ii j = k and equal to zero otherwise, then we have for a 
non-relativistic approximation, i.e., taking the limit of e — mc^ as the speed 
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of light c goes to infinity in Eq. 

c _ Tn/^S _ cg°3pj+y/g00m'2c^+(cg0ipjf-c^g00g3''pjp^ ^ 2 



5-=^ mc 



c,o.,,..-V^™c^(i-.i!;a^-4^) ^ (48) 



mc 



That is, the non-relativist ic approximation gives the energy as follows; 

e — mc^ ^ --^5^^PjPk + m(j) + (jy'pj. (49) 
2m 

It is straightforward to identify the three energy terms on the right hand 
side of Eq. fH^ as the kinetic energy, a COW potential energy, and a Sagnac 
potential energy, respectively. 



3.3 Probability Current and Classical Trajectories 

Recall that for Whitham's method, we set the wave function \1/ = xe*^ where 
X = X {X) is also a wave function and 9 = 9 [X) is a real phase at each point 
X in the spacetime. Then the probability current and scalar field J are 
given by: 

J = c^+A^ = cx~^Ax. (50) 

From Eq. (HUj) . x is a solution of the equation ■y (u) x = where u^^ is a 
unit four- vector field satisfying m^m^ = 1 and where 7(m) = u^'-f^. From the 
anticommutation relation of the Dirac gamma matrices in Eq. ([3]), we have: 

= g'^u, = I {YY + 7 V) = ^ [7^(^) + 7(^)7l • (51) 
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Moreover, using again the definition 7(m) = u^'y'^, it follows easily from the 
properties of the hermitizing matrix A [Eq. (jlj)] that x^^7('w) = [7(m)x]''' A. 
We get thus from Eqs. (HOD. (|50|) and (ISTj): 

= J^. (52) 

That is, J'^ = Ju^. Using Eqs. and together with this result, Eqs. 
(Iin])-(112]) can be written as: 

g^.u^^u" = 1, (53) 

= g^iuu", (54) 

V^(Jm'^) = 0, (55) 

g 

V^n^ - VuU^ = -Ff^i,, (56) 

where F^^/ = V/iV^ — V,yV^ is the electromagnetic field tensor. Multiply by 
u'^ and contract the index u on both sides of Eq. (l56l) . Then, using Eq. (l53il 
to set (V^Ui/) = 0, and finally raising the index we get: 

(n-V,)n'^ = ^F'^X- (57) 

Consider the integral curves x'^ (s) of the four- vector velocity field u^^ (X). 
Since the four- vector velocity field consists of unit vectors by Eq. (!53|) . the 
integral curves x'^ (s) are parameterized by arc-length s. That is, from Eqs. 
(I53D and §7^: 

ds ' 



du^ 
ds 



mc 



Note that Eqs. (l58i) are precisely the classical relativistic equations of a parti- 
cle of mass m and electric charge e in a gravitational field g^j^y in the presence 
of an electromagnetic field F^i,. Thus, the integral curves of the four- vector 
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velocity field u^{X), describing the motion of wave packets, coincide with 
the trajectories of classical relativistic particles. Note that in the absence of 
the electromagnetic field F^j., the classical trajectories f l58|) are geodesies of 
the spacetime. 

Finally, we note from Eq. fl55|) that along the integral curves of {X) 
we have: 

S = -J^.u". (59) 

This completes the proof of Theorem 2. 

Note that the wave packet equations (l29l) describe a certain congruence 
of classical trajectories ( 130|) together with a scalar density J that on each 
classical trajectory in the congruence satisfies Eq. ( |3T]) . This congruence 
satisfies certain initial conditions on a spatial submanifold M discussed in 
the next subsection. 

3.4 Mathematical Structure of the Wave Packet Equa- 
tions 

Eq. (!53l) - (l56l) can be written as follows: 



d 





= 1, 




= 9t,u'^ 




= 0, 


duo duj 


e 


dx^ dx^ 


mc 


duj duk 


e 


dx^ dx^ 


mc 



(60) 
(61) 

(62) 

^Fo,- (j = 1,2,3), (63) 



jk 



(j. A; = 1,2, 3). (64) 



Eqs. ( I60ll and ( 16T1) allow us to solve algebraically for mq, and Uj in terms 
of where j = 1, 2, 3. We will show below that Eq. (!64|) gives merely a set 
of initial conditions. Thus, we are left with only four real quasi-linear partial 
differential equations contained in Eqs. (162|) and ( 163|) . for the four real fields 
and J. Indeed, from Eq. (l36l) . we have that (!63i) and (IMl) are equivalent 
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to: 



' ^^^^ 

It follows from Eq. fj65|) that Eq. (166|) remains true for all time if and only 
if it is true at an initial time. That is because from Eq. fl^S]) . we derive: 



dx^ ydx'^j dxWx'^ dx^ \ dx^ 



(67) 



It follows that Eq. flM|l gives a set of initial conditions. Thus, provided 
that the initial conditions ( 164|) are satisfied at any initial time — i.e., on 
a spatial submanifold M — the wave packet equations (I60|) - (l64|) give rise 
to well-defined solutions. As previously stated these solutions comprise a 
congruence of classical trajectories together with a scalar density. 



3.5 Gordon Decomposition for Dirac Equations in Canon- 
ical Form 

In the WKB approximation of the DFW equation in a curved spacetime, 
classical trajectories are derived from the Gordon decomposition of the prob- 
ability current = + into a convection current and a spin current 

[3]. In this subsection, we will prove the existence of the Gordon decom- 
position for all normal Dirac equations, noting that both DFW and canonical 
equations are normal. We will further show that in the Whitham approxima- 
tion, the spin current Jf vanishes, which explains the negligible effect of spin 
on wave packet solutions. More specifically, we will show that wave packet 
solutions of the form \1/ = x^*^, where x is slowly changing compared to a 
rapidly changing phase 9, can only exist if the spin current is negligible. 
As discussed above, this definition of wave packet is independent of the size 
of Planck's constant h. Indeed, no assumption will be made in this subsec- 
tion regarding the size of Planck's constant h or the speed of light c, both of 
which we will set equal to one, h = c= 1. 

For a normal Dirac equation, define the probability current J^, the spin 
current J^, and the convection current as follows: 
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2m 

Jl! = ^g^^''[^+A{D,^)-{D,yi/)+A^], (68) 

where cr'^'' = | {'j'^'j'^ — I'^l^) are the Dirac spin matrices. Substituting the 
normal Dirac equation dHj) written in the following form: 

m = —7^/}^^ (69) 
m 

into the formula for J'^ in Eq. fl^Sjl . using the gamma matrix formula 7'^7'^ = 
g'^" — ia^'^, and noting the asymmetry a'^'^ = —a^^, we get from Eqs. ( 168|) 
and (IMD: 



^ ^+^7'^^ = ^ [^+A7^7'^(D^*) - (D^$)+A7V^] 

= ^ (^^'^ - za^'^) (Z^^^) - (D^^)+A (^^'^ + za'^'^) ^] 

= J^ + J^- (70) 

That is, = Jf! + Jj^, which proves the Gordon decomposition for normal 
Dirac equations. 

Note that the probability current = + Jf is covariantly conserved. 
Consequently, the probability density current \/—gJ^ is conserved. In gen- 
eral, the currents and are not separately covariantly conserved, unless 
the coefficient fields (7^^, A) are covariantly constant. However, it is worthy 
to note that for Dirac equations transformed into canonical form, we may re- 
place the covariant derivatives D^, including electromagnetic field potentials 

as in Eq. (123 p . with the Levi-Civita covariant derivatives + ieV^. |fl 
Then the spin current and the convection current become: 

= — r*+Aa^"(9,^) + (d^^HyAa^""^] , 
2m 

Jc = ^g'"'[^^A{d, + zeV,)^-{{d, + ieV,)^)+A^]. (71) 
2m 



* Note that the probabihty current = '^^A-y'^'^ is invariant under local similarity 
transformations S defined in Eq. ([5]), so that when transforming a Dirac equation into 
canonical form, it is only and which change their form. 
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Note that the convection current in Eq. (171 p closely resembles the spin 
zero current of the Klein-Gordon equation in the presence of both gravi- 
tational and electromagnetic external fields [35], whereas the spin motion 
resides in the spin current |3]. Recall that we are using Whitham's ap- 
proximation, for which we set = xe*^ and neglect dfj_x with respect to 
{d^6)x [see before Eq. fl7r|) ]. Using the definitions in Eq. fl25]l . this gives us 

((9^ + ieK)* ~ i{dyO + eK)* = -imu^m . (72) 

Then we use the definitions in Eq. (125]) together with Eqs. (]7T]) and (172]) to 
obtain the following formulas for the spin current and the convection current: 

^ [^+Aa^"'u^^ - ^+Aa'"^u^^] = 0, 

^ ^g'^'u^ + ^+A^] = Juf". (73) 

Thus, in the Whitham approximation, the spin current vanishes and the 
convection current Jjf = Ju^. Therefore, the probability current equals 
the convection current J^. 

4 The Classical-Quantum Correspondence 

We proved in Theorem 2 that the solutions of the Whitham approximation to 
the Dirac equation, Eq. (129]) . consist of a four- velocity vector field whose 
integral curves are classical trajectories, and a scalar field J representing 
a conserved particle density — see Eqs. (130]) and (|3T]) . From Whitham's 
approximation of the Dirac equation, we have derived the dispersion relation 
(IH]) . the motion of wave packets along classical trajectories (15U]) . conservation 
of the probability current (13T]) . as well as the generalized de Broglie relations 
= hK^ in Eqs. (|35]) and (|36]). 

We will conclude this paper by summarizing results from a previous anal- 
ysis of the "classical-quantum correspondence" based on the dispersion re- 
lation (I44p alone, which can be applied to the Dirac equation [12]. We will 
interpret the fact that the integral curves of the four-velocity field are 
classical trajectories as a "geometrical optics limit" of the Dirac equation. 
However, as we have seen, unlike the WKB limit, this "geometrical optics 
limit" is one that places no restriction on Planck's constant. Whitham's 
Lagrangian method, applied to the Dirac Lagrangian (122]) transformed into 
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equivalent canonical form fl26p . has also many advantages over using the 
dispersion relation (jUj) alone as the starting point for a wave packet approx- 
imation IHj. Whitham's method preserves the conservation laws inherent in 
the starting Lagrangian fl26|) . and in particular, the Whitham wave packet 
equations conserve the probability current = Ju^. 

To any linear partial differential equation for scalar wave functions \E', of 
the form: 

a (X) VI/ + V a' {X) ^ ^""^^ = (74) 



n=l 



(summing over coordinate indices /i^ = 0, 1, 2, 3 for r = 1, 2, 3, ■ ■ ■ , n and over 
the index n = l,2,3,---,ci) where the coefficient fields a (X) and (X) 
depend on the spacetime point X— one may associate its dispersion polyno- 
mial IIx That is to say, a polynomial function of covector fields Ky, at 
each fixed spacetime point X is given by: 

d 

[K) = a {X) + 5^^V^ ^" {X) K,, K,^. (75) 

n=l 

The dispersion relation is thereby obtained from the polynomial equation 
Tlx {K) = at each fixed spacetime point X by solving for the time com- 
ponent Kq. Applications of this one-to-one correspondence are discussed in 

Ref. ng. 

This applies also if the wave functions \1/ have m components and the 
coefficient fields a {X) and a ^ " (X) are m x m matrices, as is the case 
for the Dirac equation [19], [3H]. Note that in the matrix case, the dispersion 
relation is obtained from the scalar polynomial equation det IIx (K) = 0. 
Consider a dispersion polynomial (!75|) where a {X) and a''^ {X) are m x 
m matrices. By solving det Ux {K) = for the component Kq we get a 
dispersion relation: 



w = iy(k,x,t) (76) 

expressing the angular frequency u = —cKq as a function of the spatial wave 
covector k = {Ki, K2, K^), the spatial coordinates x = (x^, x^, x^), and time 
t, together with the auxiliary equations [1], [T^ : 
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dKi du 

H = 

dt dx' 

= 0, (77) 



OK, OK, 



dx dx^ 

for j,k = 1,2,3. In general, none or multiple such dispersion relations ( !76|) 
can be derived as distinct real roots of the polynomial equation det IIx (w, k) = 
when solving for u. Assuming at least one real root, let us choose one of 
them to be = (k, x,t). Then from Eqs. f l76p and fl77p one derives the 
Hamiltonian system as in Ref. [12], Sect. 2.2: 

dKj __dW dx^ _ dW 

Noting in Eq. ( !76|) that W = —cKq, and recalling the generalized de Broglie 
relations = —Pfj_ = hK^ derived from wave packet motion in Eqs. (1351) 
and we see that Eq. (|7S]) leads us to define the Hamiltonian H = hW 
and the momentum variables Pj = hKj for j = 1, 2, 3, whereby a system of 
classical point particle trajectories emerges as follows: Indeed, solving the 
dispersion equation (jHj) for the energy e = H {p, x,t) where p = (Pi, P2, P3) 
is equivalent to solving it for the angular frequency u = W (k, x, t) as in Eq. 
(j75|) . Then Eq. (1751) is equivalent to: 

dPj _ dH dx^ _dH 

IT ~ ~dx^' ~~di ~ Wj' ^ 

One can show that the trajectories associated with the Hamiltonian that 
is, the solution trajectories of the Hamiltonian equations ( 179|) . are identical 
to the solution trajectories of the Euler-Lagrange equations deduced from 
the well known Lagrangian for classical point particles in a background of 
electromagnetic and gravitational fields, which is given by [39], [IQ]: 



, dx^^ dx^ e^^dx'^ 

where ^ is an arbitrary parameter for the classical trajectory {^). To show 
this one first applies an inverse Legendre transformation!^ to the Hamiltonian 



^ The Legendre transformation is its own inverse |40] . pages 563-565. Thus, an inverse 
Legendre transformation is also a Legendre transformation. 
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H = .f/" (p, x,t) to obtain a traditional Lagrangian L (x, and then, as 

in Ref. [ID], pages 267-271, one generalizes the trajectory parameter to be 
an arbitrary parameter ^, instead of the coordinate time t. It is straightfor- 
ward to check that the equations for the classical trajectories fl30l) are the 
Euler— Lagrange equations for the Lagrangian fISOl) . Thus, the dispersion re- 
lation ( 176|) and the auxiliary equations ( 1771) give rise to the classical point 
particle trajectory equations ( l30l) . as is also the case for the integral curves 
of the Whitham equations (!29|) . 

Similar to the previous analysis of the classical-quantum correspondence 
([12], Sect. 2.3), it is in the "geometrical optics limit" that the solutions of 
each Dirac equation transformed into equivalent canonical form satisfy the 
dispersion equation (jUj). Indeed, Whitham's approximation: d^x (9^9) x, 
which we applied in Section [3] to the Dirac Lagrangian (l26l) . is one way of 
defining precisely this limit. However, the classical Hamiltonian equations 
fl79|) . which are based solely on the dispersion relation, give only part of the 
Whitham equations ( !29l) . In addition to providing equations ( l30l) equivalent 
to the Hamiltonian equations (179|) . the Whitham equations ( l29l) preserve the 
symmetries of the Dirac Lagrangian ( l26l) . and provide for the conservation 
of the probability current = Ju'^, which is a property inherited from the 
exact Dirac equations in a curved spacetime. 

A Appendix: Proof of Theorem 1 

To prove Theorem 1 of Section [2] we will use the following corollary of a deep 
theorem of linear hyperbolic partial differential equations: 

THEOREM 0. Let Mi, M2, ■ ■ ■ , M„ be complex dxd matrix functions that 
depend smoothly on n-|-l independent real variables t,xi,X2, ■ ■ ■ , x„ in a slab 
-T < t < T, X G R", denoted as I x R". Suppose that Mq, Mi, M2, ■ ■ ■ , M„ 
are Hermitian matrices, and furthermore assume that Mq is positive definite. 
Let F = F [S) be a homogeneous linear function of complex dxd matrices S, 
as well as having explicit dependence on t,Xi,X2, ■ ■ ■ , Xn- Then the complex 
linear hyperbolic system: 

^of (81) 

has a smooth dxd complex matrix valued solution 3:1 x R° — t- M {C,d) 
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which equals the identity matrix at t = 0, as its prescribed smooth initial 
data. 

Proof: Theorem is Corollary 3 in Ref. [20], Sect. 6, which is based on a 
theorem of Lax [33] (see also Ref. [31] ) . 

Proof of Theorem 1 of Section 2. 

First note from Eq. (|4]) that the matrices = A'-/'^ are Hermitian matrices. 
Then, note that the conditions stated in Theorem 1 for the metric compo- 
nents g^i, imply that = A--/^ is a positive definite matrix by Theorem 6 
of Ref [22], Appendix B. By Theorem 3 of Ref [21], Sect. 3.4, Eq. (54), 
a local similarity transformation T of the second kind, takes a generalized 
Dirac equation of the form ffTTl) into a Dirac equation of the normal form 
f[T^ . if and only if T satisfies the following partial differential equation: 

B^D^T = -]^{D^B^')T. (82) 

If we can solve Eq. (1821) for such a local similarity transformation T, then 
the transformed coefficient fields will be as in Eqs. ^ and ([8]): 

A = T+AT, (83) 

f = r 

Now from Eq. f|T8|) . a local similarity transformation S of the first kind, takes 
a Dirac equation of the normal form f|T^ into a normal QRD-0 equation f|T6|) . 
if and only if S satisfies the following partial differential equation: 

B^d^S = -B^f^S, (84) 

where B'^ = Aj^. The matrices i?^ and i?^ are Hermitian, and moreover, 
the matrices B^ and are positive definite. Thus, the two systems (1521) and 
flM|l have same form as the complex linear hyperbolic system in Eq. flHTl) . 

Let x^U — )■ Rx R^, mapping X — (t, x), be the coordinate chart that we 
assume to be defined on U. Then, consider the projection map 7r:R x R^ — >■ 
R taking (t, x) — )■ t. Let Xq G U and let to = n o x (-^o) ^ R- Let 
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M = (tt o x) ^ (to)- Note that G M C U. It will suffice to prove that 
there exist nonsingular solutions T and S of the systems (1821) and that 
are both defined in a common open neighborhood W of Xq G U. 

By Theorem 0, the Cauchy problem for fl521) with the smooth initial data 
^|m = I4 has a smooth solution T in an open neighborhood W* of Xq. 
Denote by W the open subset of W* in which T is a nonsingular matrix so 
that exists. Note that Xq e W since Xq G M and T|m = l4- Thus, W 
is an open neighborhood of Xq such that the local similarity transformation 
T (and its inverse T~^) is well defined on W, and hence, from Eq. (155]) . the 
complex linear hyperbolic system fl84p is well defined on W. 

By Theorem 0, the Cauchy problem for fl8^ with the smooth initial data 
'S'Impiw = I4 has a smooth solution S in an open neighborhood W* C W of 
Xq. Denote by W the open subset of^W* in which S" is a nonsingular matrix 
so that exists. Note that G W since Xq G M W and S'Imhw = l4- 
Thus, W C W* C W is an open neighborhood of Xq such that both the 
local similarity transformation 5" and the local similarity transformation T 
(and their inverses 5*"^ and T~^) are well defined on W, and therefore the 
local similarity transformation T o S (and its inverse 5*"^ o T~^) is also well 
defined on W. Q. E. D. 



References 

[1] D. Wick, The Infamous Boundary: Seven Decades of Controversy in Quantum 
Physics, (Birkhauser, 1995), pages 24 and 50—53. 

[2] M. P. Silverman, More than One Mystery: Explorations in Quantum Interfer- 
ence, (Springer-Verlag, 1995), pages 1—8. 

[3] J. Audretsch, "Trajectories and spin motion of massive spin-1/2 particles in 
gravitational fields", J. Phys. A: Math. Gen. 14, 411-422 (1981). 

[4] L. E. Ballentine, Quantum Mechanics: A Modern Development (World Scien- 
tific, 2001), pages 230 and 388-390. 

[5] D. Park, Classical Dynamics and its Quantum Analogues (Springer Verlag, 
1990), pages 49-51. 



29 



[6] P. Schliiter, K. H. Wietschorke, and W. Greiner, "The Dirac equation in or- 
thogonal coordinate systems: I . The local representation", J. Phys. A: Math. 
Gen. 16, 1999-2016 (1983). 

[7] P. Schliiter, K. H. Wietschorke, and W. Greiner, "The Dirac equation in or- 
thogonal coordinate systems: II . The two center Dirac equation", J. Phys. A: 
Math. Gen. 16, 2017-2034 (1983). 

[8] G. C. McVittie, "Dirac Equation in General Relativity", Royal Astronomical 
Society (1932), pages 868-877. 

[9] G. B. Whitham, Linear and Nonlinear Waves, (Wiley-Interscience, 1995), 
pages 363-402. 

[10] F. W. Hehl and W. T. Ni, "Inertial effects of a Dirac particle", Phys. Rev. D 
42, 2045-2048 (1990). 

[11] J. Kessler, Polarized Electrons, (Springer-Verlag, 1976), pages 2—6. 

[12] M. Arminjon, "On the relation Hamiltonian — wave equation, and on non- 
spreading solutions of Schrodinger's equation", II Nuovo Gimento 114B, 
71-86 (1999). 

[13] W. Pauli, "Contributions mathematiques a la theorie des matrices de Dirac", 
Ann. Inst. Henri Poincare 6, 109—136 (1936). 

[14] W. Kofink, "Zur Mathematik der Diracmatrizen: Die Bargmannsche Hermi- 
tisierungsmatrix A and die Paulische Transpositionsmatrix B", Math. Z. 51, 
702-711 (1949). 

[15] H. Weyl, "Elektron und Gravitation", Z. Phys. 56, 330-352 (1929). 

[16] V. A. Fock, "Geometrisierung der Diracschen Theorie des Elektrons", Z . 
Phys. 57, 261-277 (1929). 

[17] D. R. Brill and J. A. Wheeler, "Interaction of Neutrinos and Gravitational 
Fields", Revs. Mod. Phys. 29, 465-479 (1957). 

[18] T. C. Chapman and D. J. Leiter, "On the generally covariant Dirac equation", 
Am. J. Phys. 44, 858-862 (1976). 

[19] M. Arminjon, "Dirac-type equations in a gravitational field, with vector wave 
function", Found. Phys. 38, 1020-1045 (2008). [arXiv:gr-qc/0702048v4] 



30 



[20] M. Arminjon and F. Reifler, "Four-vector vs. four-scalar representation of the 
Dirac wave function", Int. J. Geom. Methods Mod. Phys. 9, No. 4, 1250026 
(2012) [23 pages]. [arXiv:1012.2327v2 (gr-qc)] 

[21] M. Arminjon and F. Reifler, "Basic quantum mechanics for three Dirac 
equations in a curved spacetime", Braz. J. Phys. 40, 242—255 (2010). 
[arXiv:0807.0570v3 (gr-qc)] 

[22] M. Arminjon and F. Reifler, "Dirac Equation: Representation Indepen- 
dence and Tensor Transformation", Braz. J. Phys. 38, 248-258 (2008). 
[arXiv:0707.1829v3 (quant-ph)] 

[23] K. Godel, "An Example of a New Type of Cosmological Solutions of Einstein's 
Field Equations of Gravitation", Revs, of Mod. Phys. 21, 447-450 (1949). 

[24] M. Arminjon and F. Reifler, "General reference frames and their associ- 
ated space manifolds". Int. J. Geom. Methods Mod. Phys. 8, 155—165 (2011). 
[arXiv:1003.3521v2 (gr-qc)] 

[25] M. Arminjon and F. Reifler, "A non-uniqueness problem of the Dirac 
theory in a curved spacetime", Annalen der Physik 523, 531-551 (2011). 
[arXiv:0905.3686 (gr-qc)] 

[26] M. V. Gorbatenko and V. P. Neznamov, "Solution of the problem of unique- 
ness and Hermiticity of Hamiltonians for Dirac particles in gravitational fields" , 
Phys. Rev. D 82, 104056 (2010). [arXiv: 1007.4631 (gr-qc)] 

[27] M. V. Gorbatenko and V. P. Neznamov, "Uniqueness and self-conjugacy of 
Dirac Hamiltonians in arbitrary gravitational fields", Phys. Rev. D 83, 105002 
(2011). [arXiv:1102.4067vl (gr-qc)] 

[28] A. O. Barut and I. H. Duru, "Exact solutions of the Dirac equation in spatially 
flat Robertson- Walker space-times", Phys. Rev. D 36, 3705-3711 (1987). 

[29] G. V. Shishkin and V. M. Villalba, "Neutrino in the presence of gravitational 
flelds: Exact solutions", J. Math. Phys. 33, 4037-4045 (1992). [arXiv:hep- 
th/9307061] 

[30] V. M. Villalba, "Dirac spinor in a nonstationary Godel type cosmological 
universe". Mod. Phys. Letters A 8, 3011-3015 (1993). [arXiv:gr-qc/9309019] 

[31] E. Montaldi and A. Zecca, "Neutrino Wave Equation in the Robertson- Walker 
Geometry", Int. J. Theor. Phys. 33, 1053-1062 (1994). 



31 



[32] R. Portugal, "Exact solutions of the Dirac equations in an anisotropic cosmo- 
logical background", J. Math. Phys. 36, 4296-4300 (1995). 

[33] P. D. Lax, Hyperbolic Partial Differential Equations, (Am. Math. Soc, 2006), 
page 61. 

[34] R. Geroch, "Partial Differential Equations of Physics", arXiv:gr- 
qc/9602055vl (1996). 

[35] L. Ryder, Introduction to General Relativity, (Cambridge University Press, 
2009), pages 397, 416-418. 

[36] B. Mashlioon, "Neutron Interferometry in a Rotating Frame of Reference", 
Phys. Rev. Letters 61, 2639-2642 (1988). 

[37] L. Ryder, "Spin-rotation coupling and Fermi- Walker transport", Gen. Rel. 
Gravit. 40, 1111-1115 (2008). 

[38] M. Arminjon, "Dirac Equation from the Hamiltonian and the Case with a 
Gravitational Field", Found, of Phys. Letters 19, 225-247 (2006). [arXiv:gr- 
qc/0512046v2] 

[39] B. L. Moiseiwitsch, Variational Principles, (Dover, 2004), pages 45—50. 

[40] O. D. Johns, Analytic Mechanics for Relativity and Quantum Mechanics, (Ox- 
ford University Press, 2005), pages 267—271 and 563—565. 



32 



